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I. INTRODUCTION 



It has been realized in the recent years that a palette of 
nanocharacterization tools has to be employed in order to 
elucidate the origin of the surprising high-temperature ferro- 
magnetism detected in a number of magnetically doped semi- 
conductors and oxides. The application of state-of-the 
art analytic methods appears to reveal that the distribution 
of magnetic cations is generally highly non-uniform and that 
there is an ultimate relationship between the ion arrangement 
and the magnetic response of the system. In particular, ran- 
domly distributed localized spins account for the paramag- 
netic component of the magnetization, whereas regions with 
a high local density of magnetic cations are presumably re- 
sponsible for ferromagnetic features. It has been found that 
the aggregation of magnetic ions - driven by a substantial 
contribution of open d shells to the bonding - leads either 
to crystallographic phase separation, i.e., to the precipitation 
of a magnetic compound or an elemental ferromagnet, or to 
the chemical phase separation into alternating regions with 
higher and lower concentration of magnetic cations, occur- 
ring without distorting the crystallographic structure. In the 
literature on diluted magnetic semiconductors (DMS), this 
chemical phase separation is commonly referred to as spin- 
odal decomposition,''^ independently of its microscopic ori- 



gin. 

There are already experimental indications that the above 
scenario applies to the case of (Ga,Fe)N fabricated by means 
of metalorganic vapor phase epitaxy (MOVPE),^ '* '' where - 
depending on the growth parameters and on the Fe concen- 
tration - Fe is either randomly distributed over cation sites 
or non-uniformly incorporated, showing spinodal decomposi- 
tion or the presence of precipitates that have been identified 
as hexagonal e— FesN, according to synchrotron radiation x- 
ray diffraction (XRD) and high-resolution transmission elec- 
tron microscopy (HRTEM).^ '^ Moreover, by employing these 
nanocale characterization methods it has been found that the 
aggregation of Fe cations is largely hampered by co-doping 
with Si donors or Mg acceptors.^ 

In this work, we present results of x-ray absorption fine 
structure (XAFS) spectroscopy '° carried out on (Ga,Fe)N and 
(Ga,Fe)N:Si samples. XAFS is a well established tool in the 
study of semiconductor heterostructures and nanostructures" 
and has proven its power as a chemically sensitive local 
probe for the site identification of Mn/Fe dopants in III-V 
12,13,14,15,16 ^jjjj Mn-rich nano-columns in Ge.'^ In 
the present case it allows us to obtain information on the lo- 
cal atom arrangement in the vicinity of Fe ions as well as 
on the Fe charge state, depending on the Fe concentration, 
growth conditions, and co-doping. We then compare the ex- 
perimentally determined bond lengths to the corresponding 
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values expected from ab initio simulations. This combined 
experimental and computational effort makes it possible to 
identify the epitaxial parameters controlling the formation of 
e-FeaN precipitates in (Ga,Fe)N. However, no presence of 
bonds specific to e-Fe3N has been detected in the (Ga,Fe)N 
films co-doped with Si. This finding provides a strong cor- 
roboration of the previous theoretical suggestions'** and ex- 
perimental observations^ indicating that the alteration of the 
magnetic ions valence by shallow impurities can hinder their 
aggregation and, therefore, extend the solubility range. 

Our paper is organized as follows. In Sec. II A we present 
the relevant characteristics of the studied samples and we 
briefly recall the employed growth and doping procedures. 
The XAFS method and the analysis of both the Extended (EX- 
AFS) and Near-Edge (XANES) regions of collected spectra 
are described in Sec. II B. The employed ab initio compu- 
tation scheme are summarized in Sec. II C for the case of 
wurtzite (Ga,Fe)N and e-FesN. Section III contains an inter- 
pretation of the experimentally determined values of the bond 
lengths in the light of the theoretical results. As discussed in 
Sec. IV, this insight allows to quantify the relative concentra- 
tion of Fe in the e-FcaN precipitates depending on the growth 
conditions and co-doping with Si. The main conclusions and 
outlook of our work are summarized in Sec. V. 



II. EXPERIMENTAL 

A. Samples production 

The (Ga,Fe)N and (Ga,Fe)N:Si epilayers considered here 
are fabricated by MOVPE on c-plane sapphire substrates fol- 
lowing the growth procedures and applying the in situ and ex 
situ characterization methods we have formerly reported for 
(Ga,Fe)N.^ The total Fe concentration in the samples ranges 
from 4x10^^^ cm^'^ to 3x10^" cm^'^ for a Fe-precursor 
(Cp2Fe) flow-rate from 50 to 350 standard cubic centimeters 
per minute (sccm).^ In the Si-doped (Ga,Fe)N structures the 
Si-content is estimated to be 1 x lO'^^ cm"-^. The growth-rate 
during deposition is regulated by the Ga-precursor (TMGa) 
flow-rate and varies from 0.08 nm/s for 5 seem to 0.3 nm/s for 
12 seem of TMGa flow. 



B. Experiments and data analysis 

The XAFS measurements at the Fe-K edge (7112 eV) 
are carried out at the "GILDA" Italian collaborating research 
group beamline (BM08) at the European Synchrotron Radia- 
tion Facility in Grenoble. The monochromator is equipped 
with a pair of Si( 11 1) crystals and run in dynamical focusing 
mode.^° Harmonics rejection is achieved by using a pair of 
Pd-coated mirrors with an estimated cutoff of 18 keV. Data 
are collected in the fluorescence mode (normal geometry) us- 
ing a 13 -element hyper pure Ge detector and normalized by 
measuring the incident beam with an ion chamber filled with 
nitrogen gas. In order to minimize the effects of coherent 
scattering from the substrate, the samples are mounted on a 



vibrating sample holder and measurements are carried out 
at liquid nitrogen temperature to reduce thermal disorder. For 
each sample the integration time for each energy point and 
the number of acquired spectra are chosen in order to collect 
« 10^ counts on the final averaged spectrum. In addition, be- 
fore and after each measurement a metallic Fe reference foil 
is measured in transmission mode to check the stability of the 
energy scale and to provide an accurate calibration. In this 
way we locate at 7112.0 eV the first inflection point of the 
absorption spectrum. 

The EXAFS signal - x(fc) - is extracted from the absorp- 
tion raw data - ij,{E) - with the VIPER program^^ employing 
a smoothing spline algorithm and choosing the energy edge 
value (Eq) at the maximum derivative. The quantitative anal- 
ysis is caiTied out with the IFEFFIT-ARTEMIS programs^"*'^^ 
with the atomic models described below. Theoretical EXAFS 
signals are computed with the FEFF8.4 code^^ using muffin 
tin potentials and the Hedin-Lunqvist approximation for their 
energy-dependent part. The free fitting parameters used in 
the analysis are: (common amplitude parameter), AEq 
(the refinement of the edge position), Ri and af, the inter- 
atomic distance and the Debye-Waller factor for the z*'' atomic 
shell around Fe respectively. In addition, a linear combina- 
tion parameter {X) is fitted when two theoretical models are 
used in the same fit. The fits are carried out in the Fourier- 
transformed space (R space) in the range [1.3-3.5] A from 
file fc2 -weighted EXAFS data in the range [2.5 - 10.5] A^^ 
using Hanning windows with slope parameter dk ~ 1 and 
dR = 0.1, for the forward and backward Fourier transforms, 
respectively. In order to reduce the variables used in the fit 
and keep the theoretical models coherent for all samples, the 
fits are limited to i? = 3.5 A. In fact, the substitutional site in 
the wurtzite lattice can be reduced to two average distances of 
Fe-N and Fe-Ga and for e-FesN to the two average distances 
Fe-N and Fe-Fe. In addition, for simplicity, only one average 
Fe-N distance is reported in the text. On the other hand, for 
samples that present only the substitutional phase it is possi- 
ble to expand the model to upper distances including multiple 
scattering (MS) contributions (as pointed out, for example, in 
the Wurtzite phase analysis section of Ref. 27) but the results 
are found to be equivalent. In the case of mixed phases, the 
expansion to longer distances and MS paths, reduces consid- 
erably the overall quality of the fits due to the lack of a model 
compound for e-FesN. 

The XANES spectra are normalized using the ATHENA 
program,^^ setting the edge jump value to unity. The peaks 
appearing in the energy region before the edge are analyzed 
with the FITYK program^** to a curve consisting in an arctan- 
gent background^^ plus one or two pseudo-Voigt peaks. 



C. Theoretical methods 

The properties of a Fe atom substituting a Ga atom in 
GaN (FcGa) as well as of the e-FesN hexagonal phase 
are investigated by density functional theory (DFT) meth- 
ods by using both the local spin density - generalized gra- 
dient approximation (LSD-GGA)^^ j^^^j j^e LSD-GGAh-U 
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formalism as implemented with plane wave basis sets in 
QUANTUM-ESPRESSO,-''*'''^'-"' in order to take into account the 
strong localization of the d states of Fe, poorly described 
by LSD-GGA exchange-correlation functionals. Total en- 
ergies are calculated in a supercell approach, by using ul- 
trasoft pseudopotentials^^ planewave basis sets, the special- 
points technique for k-space integration, and the Perdew- 
Burke-Ernzerhof (PBE) exchange-correlation functional. 

In detail, for Feca the (1,1,1) k-point Monkhorst-Pack 
mesh for a 72-atoms wurtzite supercell of GaN (correspond- 
ing to 3x3x2 unit cells), a gaussian smearing of the occupation 
numbers and plane waves cutoffs of 25 Ry for wavefunctions 
and 150 Ry for densities are used. The electronic channels 
considered in the atomic pseudopotentials are 2s and 2p for N, 
3s, 3p, 3d, 4s, 4p for Fe, and 3d, 4s, 4p for Ga. One substitu- 
tional Fe atom is included in the 72-atoms supercell of GaN. 
A neutralizing background charge is imposed when dealing 
with charged states of Feca- Geometry optimizations are 
performed by fully relaxing the positions of all the supercell 
atoms by minimizing the atomic forces. The spin state of the 
system is self-consistently determined during the wavefunc- 
tion optimization. The position of electronic levels induced 
by the Fe impurity in the GaN energy gap is estimated by 
calculating the corresponding transition energy levels 
that is, the Fermi-energy values for which the charge of the 
defect changes from n to n + 1; this gives the position of 
the Fermi energy (Ep) with respect to the top of the valence 
band. These values are estimated from total energies (defect 
formation energies) as described in Refs. 38,39, where further 
details on the theoretical methods can be found. Transition en- 
ergy values have to be located with respect to the GaN energy 
gap estimated here by the e^^^ transition level relative to bulk 
GaN, thus permitting to compare defect transition levels with 
an energy gap calculated in a consistent way.'*'^ Kohn-Sham 
electronic eigenvalues at the F point and electronic density of 
states (DOS) have been also considered when discussing the 
electronic properties of the Feca impurity. 

In the case of e-FeaN, the starting point of the calculations 
is the hexagonal phase as found in literature on the nitridation 
process of Fe.'*'''*^ The space group is P6322 (No. 182) with 
the unit cell composed of 6 Fe atoms in the Wickoff site 6g 
[{x, 0, 0); (0, X, 0); {-x, -x, 0); (-x, 0, 1/2); (0, -x, 1/2); (x, 
X, 1/2)] where x = 0.333 and 2 N in site 2c [(1/3, 2/3, 1/4); 
(2/3, 1/3, 3/4)], that is, iron atoms show the motif of a slightly 
distorted hexagonal close packing (hep) structure and nitro- 
gen atoms occupy only corner-sharing octahedra (Fig. 7 in 
Ref. 42). Satisfactorily convergent results are achieved by us- 
ing the (12,12,8) k-point Monkhorst-Pack mesh, plane waves 
cutoffs of 35 Ry for wavefunctions and 140 Ry for densities. 



m. RESULTS 

A. Solubility limit as a function of Fe content 

The first series systematically studied consists of samples 
prepared each with different Fe content (Cp2Fe from 150 to 
350 seem) at a fixed growth-rate (TMGa = 5 seem). The re- 
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FIG. 1: (Color online) fc^ --weighted EXAFS signals for (Ga,Fe)N 
samples at TMGa = 5 seem as a funetion of Cp2Fe content (seem, 
labels on spectra); vertical dashed lines highlight some parts of 
the spectra as described in Sec. Ill A; kmin and kmax delimit the 
Fourier-transformed part of the spectrum used in the fit. 
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FIG. 2: (Color online) Amplitude of the FT for spectra shown in 
Fig. 1 with relative fits . Rmin and Rmax delimit the fit region; _Ri, 
R2 and Rz (with relative dashed vertical lines) indicate the main av- 
erage distances found in the fit. The R scale has no phase correction, 
so that all the peaks appear shifted by ~ 0.4 A. 



cent investigation by HRTEM shows on this series a solu- 
bility limit at 200 seem at our growth conditions. This means 
that after this limit precipitates appear mostly as e-FcsN and 
only at the surface of selected samples as a-Fe. In Fig. 1 are 
plotted the fc^ -weighted EXAFS data where qualitative differ- 
ences are clearly visible and are highlighted by the peaks A, 
B and C. In particular, three different spectra are visible: for 
150 seem, for 200 and 250 seem, for 300 and 350 scmm. This 
difference is also found in the Fourier transformed (FT) spec- 
tra shown in Fig. 2 where, in the range Rmin-Rmax, from a 
two-peaks situation (i?i and R^) at 150 seem, an intermedi- 
ate third peak {R2) appears and increases in amplitude with 
increasing Fe content. 
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TABLE I: Results of the EXAFS quantitative analysis for the first (Ga,Fe)N series at TMGa = 5 seem (A), the seeond (Ga,Fe)N series at Cp2Fe 
= 300 seem (B) and the third (Ga,Fe)N:Si series as a funetion of both TMGa and Cp2Fe (C); for elarity only average distanees are reported (see 
text). The values of Sq and AEo are stable, in the limit of error bars, at 0.9(1) and 6(3) eV, respeetively. Averaged erystallographie distanees 
are also reported from literature (D) for eomparison. The last part of the table (E), report LSD+U atomie distanees ealeulated as deseribed in 
the text. Errors on the last signifieative digit are given in braekets.**' 



Cp2Fe TMGa 


<Fe-N> 


<Fe-Fe> 


<Fe-Ga> 


X 




Ri ai 




0-2 








(seem) (seem) 


(A) (lO'^A-^) 


(A) 


(lO'^'A-^) 


(A) 




(%) 


A. (Ga,Fe)N as a function of Fe flow-rate 


150 5 


2.02(1) 6(2) 






3.18(1) 


6(2) 


100(10) 


200 5 


1.98(1) 7(1) 


2.71(3) 


19(5) 


3.19(1) 


7(1) 


70(10) 


250 5 


1.98(1) 7(1) 


2.70(3) 


22(5) 


3.19(1) 


7(1) 


70(10) 


300 5 


2.01(2) 8(1) 


2.74(3) 


21(2) 


3.21(2) 


8(1) 


50(10) 


350 5 


2.02(3) 7(2) 


2.75(4) 


21(2) 


3.22(2) 


7(2) 


40(10) 


B. (Ga,Fe)N as a funetion of Ga flow-rate 


300 5 


2.01(1) 9(1) 


2.75(2) 


28(2) 


3.23(1) 


9(1) 


50(10) 


300 8 


1.98(2) 12(4) 






3.19(1) 


9(1) 


100(10) 


300 10 


1.97(1) 8(3) 






3.19(1) 


7(1) 


100(10) 


300 12 


1.97(1) 13(5) 






3.19(1) 


8(1) 


100(10) 


C. (Ga,Fe)N:Si 


300 5 


1.99(2) 13(3) 






3.19(1) 


9(1) 


100(10) 


300 10 


1.98(2) 9(2) 






3.18(1) 


7(1) 


100(10) 


250 10 


1.99(2) 14(3) 






3.18(1) 


8(1) 


100(10) 


100 10 


2.01(3) 11(4) 






3.19(1) 


9(1) 


100(10) 


D. Crystallographie 


e-FeaN-*^ 


1.927(1) 


2.703(2) 








a-Fe'''* 




2.499/2.886 








GaN^^^ 


1.95 








3.18 




E. DFT 


e-FesN 


1.89(1) 


2.67(1) 










1.97(1) 






3.23(1) 






1.99(1) 






3.22(1) 




Feci 


2.05(1) 






3.21(1) 





In the quantitative fit, the data at 150 seem ean be repro- 
duced with a two shell model consisting in a Fe-N (i?i) and 
Fe-Ga (Rs) shells. The intermediate peak {R2) is obtained 
with a Fe-Fe shell. The numerical results are shown in Ta- 
ble lA. 



below that limit a further phase is detected and revealed by 
the peak i?2- Also in this case the quantitative analysis is con- 
ducted by reproducing data either with a substitutional model 
or with a combination of substitutional plus e-FesN phase. 
The results are reported in Table IB. 



B. Fe incorporation dependence on growth rate 

The second series studied permits to investigate how the Fe 
incorporation changes with the growth rate. In fact, by fix- 
ing the Fe content at 300 seem (i.e., well above the solubility 
limit estimated to be about 200 seem), the TMGa flow-rate is 
increased from 5 to 12 seem through the samples series. The 
EXAFS spectra for this series are shown in Fig. 3 and their 
respective FTs in Fig. 4. From a qualitative point of view we 
note that for high growth rate values (TMGa > 8 seem) the 
EXAFS signal consists in two main frequencies (relative to 
the distances i?i and R3 in the FT) whereas in samples grown 



C. Si co-doping 

Finally, the effect of co-doping with Si is investigated in 
a samples series with variable TMGa and Cp2Fe. EXAFS 
spectra as a function of both TMGa and Cp2Fe are reported 
in Fig. 5 with their relative FTs in Fig. 6. By considering the 
same part of the spectra highlighted in the previous two series, 
it is visible in this case that no differences emerge between 
spectra"^^ and all present a typical FcQa signal. The Fe full in- 
clusion in substitutional sites is confirmed by the quantitative 
analysis reported in Table IC. 
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FIG. 3: (Color online) fc^ -weighted EXAFS signals for (Ga,Fe)N 
samples at Cp2Fe = 300 seem as a function of TMGa flow-rate (seem, 
labels on spectra); vertical dashed lines highlight some parts of the 
spectra as in Fig 1; kmin and kmax delimit the Fourier-tranformed 
part of the spectrum used in the fit. 
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FIG. 5: (Color online) fc^-weighted EXAFS signals for (Ga,Fe)N:Si 
samples as a function of TMGa and Cp2Fe (seem, labels on spec- 
tra); vertical dashed lines highlight some parts of the spectra as in 
Figs. 1,3; kmin and kmax delimit the Fourier-tranformed part of the 
spectrum used in the fit. 
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FIG. 4: (Color online) Amplitude of the FT for spectra shown in 
Fig. 3 with relative fits . Rmin and Rmax delimit the fit region; 
-Ri , R2 and i?3 (with relative dashed vertical lines) indicate the main 
average distances found from the fit. The R scale has no phase cor- 
rection, so that all the peaks appear shifted by ~ 0.4 A. 



FIG. 6: (Color online) Amplitude of the FT for spectra shown in 
Fig. 5 with relative fits . Rmin and Rmax delimit the fit region; _Ri 
and Ri (with relative dashed vertical lines) indicate the main average 
distances as in Figs. 2 and 4. The R scale has no phase correction, so 
that all the peaks appear shifted by ~ 0.4 A. 



D. Fe charge state 

In order to get information on the charge state of Fe in 
the studied samples, the near-edge region of the absorption 
spectra (XANES) is analyzed. The amplitude and position 
of the peaks due to the partially forbidden Is 3d transi- 
tions appearing in the pre-edge region of the absorption co- 
efficient were widely investigated in literature as a function 
of the local symmetry (tetrahedral or octahedral) and valence 
state (Fe'^+ or Fe^+) in Fe compounds. The general find- 
ing was that tetrahedrally coordinated compounds exhibit a 
single pre-edge peak with an amplitude above 10 % of the to- 
tal edge jump and a position changing with the valence state 



from « 7112 eV in the case of Fe^+ to « 7114 eV in the case 

of Pe3+ 30 

Here we present a systematic extension of a previous 
investigation.' Quantitative results are reported in Table II and 
two representative fits are plotted in Fig. 7; undoped samples 
exhibit only a single peak of average amplitude 13 ± 1 % and 
positon 7113.9 ± 0.1 eV whereas the doped samples present 
the same peak at 7113.9±0.1 eV with a slightly lower average 
amplitude of 10 ± 1 % and a further peak at 7112.8 ± 0.1 eV 
of amplitude 3 ± 1 %. These data are explained by the pres- 
ence of Fe'^^ ions in the first series and a co-existence of Fe'^^ 
and Fe-^+ ions in the second case. The appearance of double 
peaks in the pre-edge region in case of co-existence of chemi- 
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TABLE II: Quantitative results of the XANES analysis. 



Cp2Fe TMGa 


peak I 


peak II 




center 


height 


center 


height 


(seem) (seem) 


(eV) 


(%) 


(eV) 


(%) 


(Ga,Fe)N 


300 5 


7113.8(1) 


8(1) 






300 8 


7113.9(1) 


12(1) 






300 10 


7113.9(1) 


14(1) 






300 12 


7113.9(1) 


17(1) 






(Ga,Fe)N:Si 


300 5 


7113.9(1) 


10(1) 


7112.9(1) 


2(1) 


300 10 


7113.9(1) 


11(1) 


7112.9(1) 


3(1) 


250 10 


7113.9(1) 


11(1) 


7112.7(1) 


3(1) 


100 10 


7113.8(1) 


9(1) 


7112.5(1) 


4(1) 



cal species at different valence states was already pointed out 
in literature^' '*'* and in our case evidences the partial reduction 
of the metal ions (Fe''+ Fe^+) upon Si addition. 



E. Theoretical results 

1. Substitutional Fe in GaN 

The local structure of a Feca in GaN is investigated 
through geometry optimization procedures both in the LSD 
and LSD+U frameworks. The estimated lattice constants of 
the 72-atoms supercell are a = 9.66 A, c = 10.49 A. De- 
tails of the atomic geometries produced by the (more reliable) 
LSD+U calculations are given in Table IE for the charge states 
0, +1 and -1 of Feca. which correspond to the Fe^^, Fs'q'^ 
and Fe^^ forms of the impurity, respectively. Present results 
show that the Feca-N bond distances for the -1 (+1) state are 
slightly longer (almost the same) with respect to those of the 
neutral state. Moreover, no significant differences are found 
between the Feca-Ga distances estimated for the above three 
charge states of Feca- A quite similar picture is provided by 
the LSD calculations. 

Regarding the electronic properties of FcQa, total spin DOS 
and spin DOS projected on the Fe atomic orbitals as given 
by the LSD+U calculations for the neutral state of the impu- 
rity are reported in Fig. 8 together with the calculated Fermi 
energy. These DOS have to be considered together with an 
estimate of 1.7 eV given by the Kohn-Sham electronic eigen- 
values for the GaN energy gap (Eg). Then, Fig. 8 shows that 
occupied electronic states induced by the substitutional Fe are 
mixed with GaN states at the top of the valence band (TVB), 
whereas unoccupied Fe states are resonant with the conduc- 
tion band. A more reliable estimate of the location of the 
levels induced by a Feca in the GaN energy gap is achieved 
by calculating the corresponding e"/"'^ transition level where 
the Fe impurity changes its charge state from to -1. The 
present evaluation of the e^/~^ level is 2.25 eV to be com- 
pared with a calculated Eg of 2.93 eV. These results favorably 



compare with experimental results locating the Fe acceptor 
level at 2.86 eV,"*'^ on a band gap of 3.50 eV for wz-GaN at 
K.^° It can be noted that an even better agreement with the 
experiment could be reached by assuming a linear scaling of 
fO/^i with the energy gap. In this case, indeed, a value of 
2.70 eV for the Fe level would correspond to the experimental 
energy gap. 

The present results agree with a deep acceptor character of 
FcQa- They seem also to rule out a possible donor behavior 
of FcQa- In Fig. 8 it is shown indeed that there are no oc- 
cupied Fe states in the GaN energy gap. Moreover, a value 
of —0.03 eV has been estimated for the e^/'^^ transition level 
indicating that a Fe'^+/^+ donor level of Feca would be reso- 
nant with the valence band. In agreement with this conclusion 
no such state has been found experimentally in the band gap 
of GaN;Fe.^' However, the present ab initio findings do not 
rule out a possibility that the Fe impurity in GaN gives rise 
to a charge transfer state, Fe'^+ + h, acting as a hole trap, as 
postulated by some of us.^-^^ 

The results of the present calculations permit also to asses 
estimate a change in the formation energies of FcQa in dif- 
ferent charge states (ilp) as a function of the Fermi energy 
Ep. For instance, for the magnitude of rip'^-fip, we estimate 
+0.40 eV and -1.35 eV for Ep equal to Eg/2 and Eg, respec- 
tively (in this case we use the experimental value of the energy 
gap Eg). This indicates that the formation of Fe^^ is pro- 
moted in intrinsic GaN, while Fe^^ is favored in the n-doped 
material. On the contrary, the value of fip^-fip is evaluated as 
+0. 10 eV and +1 .85 eV for the Fermi level at the top of the va- 
lence band and in the middle of the energy gap, respectively, 
that is, a donor behavior of Fe (corresponding to the formation 
of Fel^) is never favored. 



2. e-FesN 

The cell parameters have been calculated by minimizing the 
total energy of the system. Optimal values have been found 
at the zero value of linear fits of the total stress, correspond- 
ing to the minimum of the total energy. Such a procedure 
gives a = 4.64 A and c = 4.34 A, that is w 2% contraction 
of the experimental values a^xp = 4.7209(6) A and Cexp = 
4.4188(9) A; this being a typical effect of the LSD approxi- 
mation. The relaxation properly reproduces the experimental 
value of X that is a 0.009 contraction from the 1/3 ideal value 
of X for Fe in the hep structure. The resulting atomic distances 
reported in Table IE are slightly different from the crystallo- 
graphic ones"*^ (Table ID). The above structural minimization 
corresponds to a magnetization of 2.2 /i^/atom, in agreement 
with the experimental value found from neutron diffraction at 
4.2 K.^^ Finally, the total DOS shown in Fig. 9 gives a clear 
evidence of the metallic behavior of this system. These results 
are achieved by LSD calculations; the LSD+U scheme is also 
tested with different U values in the range 2-5 eV and no ap- 
preciable difference are found between the LSD and LSD+U 
results. 
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FIG. 8: (Color online) LSD+U total spin density of states (DOS) and 
DOS projected on the Fe atomic orbitals for a substitutional Fe in 
the neutral charge state. The DOS for Fe are magnified ten times for 
sake of clarity. The calculated Fermi energy (Ep) is 10.52 eV and is 
indicated by a vertical dashed line to respect to the top of the valence 
band (TVB). 
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FIG. 9: (Color online) LSD total spin density of states (DOS) for 
e-FeaN. The calculated Fermi energy (Ef) is 17.64 eV. 



IV. DISCUSSION 
A. EXAFS 

In the first series presented in this study the sample with a 
lower Fe content exhibits two coordination shells around Fe. 
These shells are identified as Fe-N end Fe-Ga contributions 
and the observed bond length values are in agreement with 
those calculated via DFT for a Fe ion in a Ga-substitutional 
site. Thus, we conclude that Fe substitutes Ga, as previously 
suggested by a qualitative analysis of EXAFS data.^"* The 
closeness of the theoretical estimates of the Fe-N distances 
reported in Table IE makes it difficult to exploit the Feca local 
geometry to distinguish between its possible, different charge 
states. However, the theoretical results support a deep accep- 
tor character of FcQa and locate the corresponding electronic 
level in the upper half of the GaN energy gap. Accordingly, 
the formation energy is 0.40 eV lower in energy than 
in intrinsic GaN, thus suggesting the presence of Feca in the 
neutral state. 

For samples grown with higher Fe content an additional 
shell appears, matching a coordination at 2.70 — 2.75 A. This 
value coincides with the experimental and calculated Fe-Fe 
second shell distance in e-FcaN. This suggests that Fe partici- 
pates in the form of e-FeaN. It is worth noting that in this com- 
pound the first shell coiTesponds to the Fe-N bonds exhibiting 
the same structural parameters (length and number of neigh- 
bors) as those specific to Ga-substitutional Fe in GaN. These 
results provide an independent confirmation of our previous 
studies on (Ga,Fe)N,^''*'' according to which Fe-rich nanocrys- 
tals were evidenced by neither HRTEM nor HRXRD and a 
paramagnetic response was observed by SQUID in samples 
with a low-Fe content, whereas samples with a high Fe con- 
tent showed the presence of Fe-rich nanocrystals and a ferro- 
magnetic behavior. 

From a reduction in the amplitude of the dominating Fe-Ga 
signal relative to the pure substitutional specimen it is possible 
to estimate the relative content of Fe in the two phases (GaN 
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and e-FeaN), as shown in Table I. From these data we infer 
that the flow-rate of 200 seem Cp2Fe corresponds to the onset 
of e-FesN precipitation at our growth conditions, though the 
spinodal decomposition, that is the aggregation of Fe cations 
without crystaflographic phase separation may begin at lower 
Fe contents. Significantly, no presence of e-FesN is detected 
in samples co-doped with Si. 

It is worth to underline that no evidence is found for the 
a-Fe phase that should be witnessed by a double Fe-Fe shell 
at 2.499 and 2.886 A.'*'* Actuafly, this phase was shown to be 
present only at the surface of samples grown under particular 
conditions'*^ and the overall Fe fraction in this phase presum- 
ably lies below the detection limit of EXAFS (about 15% in 
the present case). 

We also observe that our data do not point to the presence 
of Fe in interstitial sites, at least in samples containing no sec- 
ondary phases, for which the EXAFS signal is completely re- 
produced by a simple substitutional model. In the case of sam- 
ples with precipitates, their contribution could a priori mask a 
signal coming from interstitials. This is in contrast to the case 
of (Ga,Mn)As where evidence for Mn interstitials has been 
found in EXAFS studies, i^"' 

Turning to samples deposited at different growth rates we 
note that if the Fe content exceeds the previously deter- 
mined solubility limit, we find again the presence of e-FesN 
nanocrystals in the film grown at the lowest rate (TMGa = 
5 seem). For faster growth rates, however, the EXAFS spectra 
correspond exclusively to substitutional Feca in GaN, empha- 
sizing the effectiveness of the fast growth in suppressing the 
formation of segregated phases. This observation is in qual- 
itative agreement with the previous study,^^ where HRTEM 
carried out for samples grown at low and high rates (although 
at a noticeably higher Fe content) revealed the formation of 
e-FeaN only for slowly grown specimens. 

B. XANES 

The XANES spectral region is characterized by the pres- 
ence of a single pre-edge peak at about 71 14 eV in all samples 
grown without Si. In agreement with the literature data-'" -"''*^ 
this reveals that Fe assumes the Fe'^^ charge state in a tetra- 
hedral environment, and this result corroborates the conclu- 
sions derived from previous electron paramagnetic resonance 
studies^'^' and from the present DFT calculation for (Ga,Fe)N. 
Co-doping with Si results in the appearance of an additional 
peak in the pre-edge region at about 7112.8 eV, which is at- 
tributed to Fe'^^ ions in a tetrahedral environment, as already 
pointed out in our prior work.^ The fraction of the total Fe 
content in this particular valence state is below 20 %, the ma- 
jority of Fe impurities still remaining in the Fe'^+ configura- 
tion. This explains why we do not observe in the EXAFS 
data (and, in particular, in the values of the Fe-N bond length) 
any significant deviation when comparing the samples con- 
taining no Si impurities with the Si-doped ones. However, the 
co-deposition of Si hampers the aggregation of Fe and, thus. 



shifts the solubility limit to higher Fe concentrations. 

The above interpretation of the XANES findings is consis- 
tent with the ab initio results, implying that the Fe^+/Fe'^+ 
state, i.e. the e"/~i level, resides in the GaN gap 2.25 eV 
above the top of the valence band. Thus, this state should be 
occupied by electrons provided by shallow donors, such as Si. 
This conclusion is further supported by the computed forma- 
tion energy that results to be for the Fe''+ configuration by 
1.35 eV higher than for the Fe^+ case. 



V. CONCLUSIONS AND OUTLOOK 

The present EXAFS results, together with previous 
HRTEM, electron dispersive spectroscopy, and synchrotron 
XRD,^ show how the Fe incorporation can be efficiently con- 
trolled by Fe flow, growth rate, and co-doping with Si. In par- 
ticular, the aggregation of Fe cations can be minimized by in- 
creasing the growth rate and by co-doping with Si, shifting the 
solubility limit towards higher Fe content at given growth con- 
ditions. While there is agreement between the present DFT 
computations and both the previously and here observed 
Fe charge state, including its evolution with the Si concen- 
tration, a detailed understanding of how co-doping and other 
growth parameters affect the aggregation of Fe ions awaits for 
a comprehensive theoretical treatment. 

A strict correlation between magnetic properties and the 
presence of Fe-rich nanocrystals^'^ strongly suggests that the 
surprisingly robust ferromagnetism of (Ga,Fe)N as well as 
of other diluted magnetic semiconductors and oxides de- 
posited under specific growth conditions results from the self- 
organized assembly of magnetic nanocrystals, whose high 
blocking temperatures account for the survival of ferromag- 
netic features at high temperatures. This emerging insight 
provides a promising starting point for exploiting a number 
of expected functionalities of these nanocomposite semicon- 
ductor/ferromagnetic metal systems. 
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